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FOREWORD 

The  axisymmetric  perturbations  (sausage  and  hollowing  modes)  of  an  intense 
relativistic  self-pinched  electron  beam  propagating  in  a  resistive  plasma  back¬ 
ground  are  studied,  especially  for  a  beam  with  rounded  radial  density  profile. 

The  Bennett  profiles  are  assumed  for  both  the  equilibrium  beam  current 

2  2-2 

Jb(r)  =  Jb(°)  (1  +  r  /R^)  and  plasma  return  current  Jp(r)  *-f  J^CO) 

2  2-2 

(1  +  r  /Rp)  ,  where  R^  and  R^  are  the  characteristic  radii  of  the  beam  and 
plasma  return  currents  respectively.  It  is  further  assumed  that  the  electric 
conductivity  o(r)  of  the  plasma  channel  is  proportional  to  the  return  current. 
For  a  paraxial  electron  beam  with  complete  space-charge  neutralization  by  the 
ambient  plasma,  the  axisymmetric  modes  can  be  destabilized  by  the  phase  lag 
between  the  magnetic  field  and  beam  current,  even  without  the  plasma  return 
current . 

The  plasma  return  current  significantly  modifies  the  growth  rate  of  the 

2  2 

instability  such  that  the  ratio  of  plasma  to  beam  current  (-1^/1^  ■  f  Rp/R^) 
largely  determines  the  stability  character  of  the  beam.  Furthermore,  for  the 
same  fractional  current  neutralization  f,  the  modes  are  highly  unstable  for 
smaller  ratio  of  plasma  to  beam  radius  Rp/R^.  As  comparing  to  the  resistive 
hose  instability,  the  growth  rates  for  hollowing  mode  could  be  larger  than  those 
of  hose  mode,  while  the  sausage  mode  is  much  stabler  than  the  hose  mode. 
Stability  properties  are  illustrated  in  detail  for  various  system  parameters. 
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I .  Introduction 

There  has  been  considerable  interest  in  recent  years  in  the  propagation  of 

high-energy  self-pinched  relativistic  electron  beam  in  the  highly  resistive  plasma 

medium.  Naturally,  the  equilibrium  and  stability  properties  of  the  electron 

beam  become  one  of  the  important  subjects  under  investigation.  The  resistive 

hose  instability  appears  to  be  the  most  serious  magnetohydrodynamic  instability 

for  a  beam  in  a  resistive  plasma  channel  and  has  been  studied  extensively  by 
1-3 

many  authors.  Nevertheless,  the  axisymmetrie  oscillations  could  expell  the 

4 

beam  from  its  channel  and  plague  early  beam  propagation.  It  is  the  purpose  of 

this  paper  to  examine  the  importance  of  the  axisymmetrie  modes  in  comparison  with 

the  hose  mode.  The  axisymmetrie  perturbations  of  an  intense  self-pinched 

relativistic  electron  beam  propagating  in  the  plasma  background  have  been  studied 

5  6  7  8 

previously  by  many  workers  both  analytically  ’  and  numerically.  *  Recently 

1  9 

Uhm  and  Lampe  developed  a  self-consistent  theory  '  to  derive  the  dispersion 
relation  for  a  flat-topped  radial  density  profile  electron  beam.  This  paper  is 
an  improved  version  of  the  previous  study,  including  the  influence  of  the 
rounded  beam  density  profile  on  the  axisymmetrie  modes.  It  has  been  understood 
that  the  flat-topped  electron  beam  has  the  feature  of  unbounded  linear  growth 
rate  for  axisymmetrie  oscillations  because  all  the  beam  electrons  are  assumed  to 
have  the  same  natural  oscillation  frequency  which  resonants  with  the  wave.  As 
we  know  in  the  real  world,  the  beam  contains  electrons  of  all  energies  from  the 
smallest  up  to  the  maximum  possible  energy  which  is  equal  to  the  kinetic  energy 
of  the  accelerated  electron.  The  electrons  at  different  energy  level  oscillates 
with  different  frequencies  and  radii,  thus  the  unharmonic  beam  electrons  at 
different  radial  position  exhibit  different  betatron  frequency  and  resonant 
locally  with  the  waves.  Electrons  with  an  excess  of  energy  will  move  with 
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a  greater  radius  and  electrons  with  less  energy  oscillates  around  smaller' radius. 
Therefore,  each  electron  automatically  chooses  the  appropriate  phase  necessary 
for  their  interaction  with  perturbed  waves.  This  phase  mixing  is  a  damping 
mechanism  for  oscillation  which  prevents  perturbation  from  growing  to  infinity. 

From  the  thermal  equilibrium  calculations  of  the  intense  relativistic 
electron  beam,  it  was  found  that  the  radial  density  profile  for  a  nonrotating 
electron  beam  has  the  familiar  Bennett  equilibrium  as  derived  originally  by 
Bennett. 1(1  For  the  choice  of  the  thermal  distribution  function,  the  beam 
temperature  profile  and  the  axial  velocity  profile  are  uniform  in  the  radial 
direction.1*  For  the  matter  of  practical  interest,  we  assume  the  beam  flow  is 
paraxial,  such  that  all  the  beam  particles  have  nearly  the  same  axial  velocity, 
and  the  transverse  energy  is  much  smaller  than  that  in  the  axial  direction.  For 
the  type  of  beam  we  consider  here,  the  equilibrium  beam  current  can  very  well  be 
represented  by  the  Bennett  profile 

Jb(r)  =  Jb(0)/(1  ♦  r2/R^)2  (1) 

where  is  the  characteristic  beam  radius. 

The  resistive  plasma  background  responds  the  beam  current  with  a  plasma 
return  current  which  strongly  influences  the  stability  behavior  of  the  beam. 

A  small  transverse  displacement  of  the  beam  can  separate  the  beam  from  the 
return  current  which  eventually  causes  the  instability.  For  simplicity  we 
assume  the  plasma  background  supplies  enough  ions  to  neutralize  the  beam 
electrons  completely.  Although  the  radial  profile  of  the  return  current  is 
determined  by  the  electrical  conductivity  o(r)  which  relies  heavily  on  the 
properties  of  the  ambient  plasma  and  experimental  condition,  for  the  time 
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being,  it  is  assumed  throughout  this  paper  that  the  equilibrium  plasma  return 
current  Jp(r)  and  electrical  conductivity  a(r)  have  also  the  Bennett  profile 
like 


J  (r)  =-f  Jb(0)/(1  +  r2/Rp)2  (2) 

and 

a(r)  =  a(0)/(l  +  r2/R2)2  (3) 

where  the  characteristic  beam  radius  Rp  can  be  greater  or  smaller  than  R^ 
depending  on  the  nature  of  the  background  plasma  such  as  the  ambient  pressure 
and  other  physical  properties,  f  is  the  ratio  of  beam  current  to  plasma 
current  at  the  beam  axis 


f  =  -Jp(0)/Jb(0)  .  (4) 

If  the  above  mentioned  equilibrium  current  configurations  are  used,  then  the 
effective  current  neutralization  would  be  determined  by  the  ratio  of  plasma  to 
beam  current  which  has  the  value 

f  =  -I  /I.  =  f  R^/R?  (5) 

eff  P  b  P  o 

The  influence  of  the  steady  return  current  on  the  axisymmetric  perturbations 
depends  largely  on  the  extent  of  plasma  current  profile  overlapped  with  the  beam 
current.  By  varying  and  R^,  the  effective  current  neutralization  not  only 
denotes  the  measure  of  roundness  of  the  beam  but  also  represents  the  degree  of 
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overlap  between  beam  and  plasma  current  profiles.  It  is  an  important  parameter 
used  later  to  describe  the  effect  of  the  current  profiles  on  the  instability. 

In  the  sections  which  follow,  the  equilibrium  properties  of  a  self-pinched 
electron  beam  with  a  steady  plasma  return  current  are  discussed  in  Section  II, 
the  necessary  condition  for  radial  confinement  of  the  electron  beam  is  alined 

The  standard  linearized  Vlasov-Maxwell  treatments  are  summarized  in  Sec  t>  III 
An  integro-differential  eigenvalue  equation  is  obtained  for  the  genera!  -ded 
beam  case,  making  use  of  the  energy  group  model.  The  most  important  sou._e 
term,  namely  the  perturbed  beam  current  for  the  rounded  beam  is  also  derived. 
For  the  purpose  of  comparison  to  the  previous  results,  an  integral  form  of 
an  approximate  dispersion  relation  is  obtained  for  the  general  rounded  beam. 

For  the  special  case  where  the  electrical  conductivity,  the  beam  and  plasma 
return  currents  all  have  the  Bennett  profile,  numerical  results  are  obtained 
for  various  system  parameters  such  as  f,  Rp  and  R^.  Stability  properties  for 
the  sausage  and  hollowing  modes  are  illustrated  and  discussed  in  great  detail . 
Finally,  the  growth  rates  of  the  axisymmetric  modes  are  compared  to  those 
of  the  hose  instability  and  the  conclusions  are  drawn  in  Section  IV. 
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II.  Equilibrium  state 

Let  us  consider  the  steady  state  of  an  infinitely  long  intense  relativistic 

self-pinched  electron  beam  propagating  in  a  resistive  background  plasma.  The 

beam  dynamic  properties  are  described  by  the  Vlasov-Maxwell  Theory.  In  the 

following  analysis,  the  cylindrical  polar  coordinates  (r,  9,  z)  are  used  with 

the  z-axis  along  the  axis  of  symmetry.  In  steady  state  (3/3 t  =  0)  all  the 

beam  equilibrium  quantities  are  assumed  to  be  azimuthal ly  symmetric  (3/3  9=0) 

and  infinitely  long  and  uniform  in  the  axial  direction  (3/3 z  =  0).  The  beam 

12 

with  net  current  I  is  assumed  to  satisfy  the  following  relationship 

I/I  =  v/Y  =  1/17000  By  <<  1  (6) 

cL 

where  v  is  Budker's  parameter  representing  the  number  of  electrons  per 

classical  electron  radius  length  of  the  beam,  =  Be  is  the  characteristic 

beam  electron  axial  velocity,  y  is  the  relativistic  factor  of  the  beam 

+ 

particle  and  I  is  the  Alfven-Lawson  limiting  current.  The  beam  electrons 
a  2  2  2 

motion  can  be  taken  to  be  paraxial  (P2  >>  +  p0)  where  P^  denotes  the 

three  components  of  the  canonical  angular  momentum  in  i  direction  respectively. 
The  axial  velocity  is  very  large  compared  to  the  transverse  velocity  and  is 
considered  to  be  a  constant.  In  this  regime,  the  equilibrium  properties 
of  an  electron  beam  can  be  described  by  a  velocity  distribution  function  of 
the  form 


f  =  F  (H  -  Y  ®c2)  &  (P  -  Pbl 
D  0  L 


(7) 


where  H  is  the  total  energy,  P,  is  the  axial  canonical  momentum.  The  single 
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particle  constants  of  motions  in  the  equilibrium  fields  are  given  by 


„  ,  2  4  2  2.Js 

H  =  (m  c  +  c  p  ) 


Pz  “  pz  '  e  Az  (r)/c 


where  y  =  (1  -  6  )”'i  and  8  =  v. /c  are  the  standard  relativistic  quantities, 
m  and  vfe  are  the  rest  mass  and  the  velocity  of  the  beam  electron,  c  is  the 
velocity  of  light  in  vacuum,  and  are  constants,  A^(r)  is  the  axial 
component  of  vector  potential  for  the  equilibrium  azimuthal  self-magnetic 
field.  There  is  no  externally  applied  magnetic  field  in  the  self-pinched 
beam  case  so  that  P0  =0.  As  a  result,  we  can  assume  that  f^  is  not  a 
function  of  Pa  as  shown  in  (7) .  For  the  nonrotating  beam,  the  vector 

O 

potential  A  (r)  has  only  the  axial  component  A^Cr)  which  can  be  determined 
self-consistently  from  the  z  component  of  the  Ampere's  law 


ih  (lV  =  "4ireB  [ n^r)  -  np(r)] 


where  n^Cr)  is  the  radial  number  density  distribution  of  the  electron  beam, 

np(r)  is  not  the  actual  number  density  but  related  to  the  plasma  return 

current  by  J  (r)  =  -e8c  n  (r) .  and  by  definition 
P  P 


Be  *  (9) 

where  in  deriving  Eq.  (8)  the  paraxial  beam  approximation  has  been  used.  All 
the  electrons  have  the  same  axial  velocity  Vz  =  8c  and  the  transverse  velocity 
of  the  beam  electron  is  small  compared  with  the  axial  velocity.  Eq.  (8)  can 
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be  used  to  calculate  equilibrium  azimuthal  self-magnetic  field  self-consistently 
for  any  net  current  in  the  system.  The  relationship  can  be  further  related 
to  .the  squared  betatron  frequency  as 


<*)g(r)  =  -eSB0(r)/Ybmr 


(47re2B2/Ybmr 2)  JQ  dr  r  ^(r)  -  np(r)J 


where  is  the  usual  betatron  frequency  which  characterizes  the  electron 
oscillations. 

If  the  beam  and  plasma  return  currents  have  the  assumed  Bennett  profile 
as  in  Eqs.  (1)  and  (2),  then  Eq.  (10)  gives  the  value  after  integration 


"6  M  -  °'5  “pb8b  [ 


pbBb  [l/(l*r2/R2)  -  f/ (l+r2/R2)J 


2  u 

where  wpb  =  (410^(0)6  /Ybm)  is  the  beam  plasma  frequency  at  axis.  It  is 
important  to  point  out  that  is  a  function  of  r  except  in  the  special 
case  of  a  flat-topped  beam  which  has  the  constant  value;  i.e., 


“s  ■  °'5  “pbSb  »-« 


This  is  the  most  important  distinction  between  the  rounded  and  flat-topped 
beams,  which  gives, by  no  surprise,  two  complete  different  pictures  in  the 
linear  perturbation  analysis.  In  short,  the  growth  rate  for  axisymmetric 
perturbations  goes  to  infinity  as  !1  +  u.  for  the  flat-topped  beam,  because 
the  pinch  force  is  harmonic.  All  the  beam  electrons  are  able  to  resonate 
with  the  wave  at  the  same  frequency.  On  the  other  hand,  because  the  pinch 
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force  is  anharmonic  for  the  rounded  beam,  the  electrons  have  a  wide  variety 

of  betatron  frequencies  which  could  damp  the  oscillation  by  phase  mixing 

mechanism.  It  could  have  only  radially  localized  wave-particle  resonance, 

so  that  the  growth  ?*ate  is  finite,  which  will  become  more  obvious  later. 

As  a  matter  of  fact,  if  the  squared  betatron  frequency  is  negative  at 

certain  radius  r,  then  the  equilibrium  fails,  and  the  beam  is  blown  apart. 

Therefore,  the  radial  confinement  for  the  equilibrium  state  requires  at 
2 

least  Wg(r)  >  0  for  any  radius  r  inside  the  beam.  This  assures  that  the 

magnetic  focussing  force  is  pinching  the  beam;  otherwise,  the  beam  will  be 

hollowing  out  and  eventually  will  become  an  annular  beam.  Since  we  have 

assumed  in  the  beginning  that  the  equilibrium  beam  current  has  the  Bennett 

profile,  the  minimum  requirement  for  such  equilibrium  profile  to  exist  needs 
2 

o).(r)  >  0,  which  leads  to 

P 

f  <  (1  ♦  r2/R2)/(l  +  r2/R^)  (13) 

For  the  case  of  R^  >  R^,  the  fractional  current  neutralization  f  is  bounded 
by 


f  <  1 


Rp<  Rj,,  (14) 


so  that  the  beam  current  J^r)  is  peaked  on  axis.  However,  if  R^  <  P.^ 
we  obtain  from  Eq.  (13)  that 


f  <  O^/Rp) 2 


V  V  (15) 


8 


NSWC  TR  83-14 


which  is  equivalent  to  feff  <  lm  Note  that  in  the  limit  where  R^  =R^, 
both  Eqs.  (14)  and  (15)  yield  the  same  condition,  namely  f  <1,  which  alsb 
applies  to  the  flat-topped  beam  case. 

Next  let  us  consider  an  electron  beam  with  a  "loss-cone"  distribution 
function 


fb  -  (n^r^m)  6  (H  -  Y0nc2)  5  {? z  -  Pfc) . 


(16) 


Integration  of  (16)  in  the  £  space,  according  to  definition,  leads  to  a  density 
distribution  with  flat-topped  radial  profile 


Vr) 


for  0^r<  R^ 
for  R^  <  r 


(17) 


where  n^tr)  is  the  equilibrium  beam  density  profile  and  is  defined  as  the 
effective  beam  radius  such  that  the  number  of  electrons  per  unit  length  of 
the  beam  has 


J*  00 

Nb  5  2»jf  dr  r  nb(r)  .  ,  (18) 

Note  that  the  beam  with  flat- topped  or  Bennett  profiles  carry  essentially  the 
same  amount  of  total  current,  which  is  very  important  when  we  compare  the 
growth  rates  in  both  cases  later  on. 

In  the  energy  group  model  a  beam  has  been  subdivided  into  groups  of 
particles  with  the  same  value  of  energy.  By  use  of  Eq.  (7)  the  beam  current 
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profile  can  be  expressed  by  definition  as 


J.  (r)  =  2ity.  me6c  /  dU  F  (U  ♦  y.  me' }  (19) 

b  b  J  *j(r)  b 

2 

where  U=  H  -  Y^mc  is  the  effective  transverse  energy  and  4^  (r)  =  e8Az(r), 
which  appeares  in  the  lower  limit  of  integration,  is  the  minimum  possible  value 
of  U  necessary  for  an  electron  at  r.  Thus,  the  beam  can  be  considered  as  a 
superposition  of  many  components  with  different  transverse  energy 


/oo 

dU  F  (U  +  Ybmc2)  J0(r,  U) 


where  each  component  has  the  equilibrium  current  density  profile  like 


J0(r,  U)  -{ 


-2itYbme6c 


r  <  R  (U) 
r>  R  (U) 


where  R(U)  is  the  maximum  radius  for  an  electron  with  eneTgy  U. 
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III.  Stability  Analysis 

A.  Linear  Perturbations 

In  this  section,  we  use  the  linearized  Vlasov-Maxwell  equations  to 
investigate  the  axisymmetric  perturbations  (m  =  0)  where  m  is  the  azimuthal 
harmonic  number.  We  follow  a  normal  mode  approach  by  asstuning  that  all  the 
axisymmetric  perturbed  quantities  vary  with  z  and  t  according  to 

6$(x,t)  =  6$  (r)  exp  j  (urr  +  flz/v2)J  (22) 

where  the  oscillating  angular  frequency  w  is  assumed  to  be  complex  with 
lm(uj)  >  0.  The  variables  t  and  SI  are  defined  by 

Tit-  Z/V2 

and  S2  =  u)  -  kv  . 

Z> 

Where  z  and  t  are  the  axial  position  and  time  in  the  laboratory  frame, 
k  is  the  axial  wavenumber.  Since  all  the  beam  electrons  have  nearly  the  same 
axial  velocity  for  the  paraxial  beam  we  consider  here,  it  is  more  convenient 
to  use  z  and  T  as  independent  variables  than  the  normal  use  of  z  and  t.  Thus 
SI  is  the  Doppler-shifted  eigenfrequency  seen  by  a  beam  electron,  V^x  is  a 
measure  of  the  distance  of  a  particular  segment  from  the  head  of  the  beam. 

If  each  beam  segment  is  taken  to  oscillate  at  a  fixed  real  frequency  SI  then 
the  solution  for  complex  w  can  be  solved  mathematically  from  the  dispersion 
relation.  Physically  Im(w)  represents  the  growth  rate  of  the  wave  as  we  move 
backward  from  the  head  of  the  beam. 
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IS 


After  Fourier  decomposing  the  Maxwell  equations,  it  is  straight¬ 
forward  to  show  that  for  the  case  of  long  wave  length  (k  <<  1)  and  low 
frequency  (w  «  c)  regime  the  eigenvalue  equation  for  the  perturbed  field 
6<J>  has  the  form 


1  9  9  47raui  t , 

?  s  r  if  ‘  — r  ■ 


-  T  sVr) 


(23) 


where  a  is  the  electrical  conductivity  of  the  plasma  channel  and  the  Vlasov 
perturbed  beam  current  contribution  6J^(r)  is  defined  as 


/ 


SJb(r)  =  -  e  /  vz  6fb  (r,  £)  d£ 


(24) 


where  6fb(r,£)  is  to  be  calculated  from  the  integration  of  the  linearized 
Vlasov  equation.  According  to  the  energy  group  model  in  Eq.  (20)  one  can 
express  the  perturbed  beam  current  as 


6Jb(r) 


■/ 


dU  F(U  Yb  me  )  5J0(r,U). 


(25) 


By  use  of  Eq.  (19)  one  can  derive  the  following  relationship 


Fd$0(R)/ dR  *  (l/27TYbmeSc)  9Jb(R)/9R 


(26) 


Substituting  Eq.  (26)  into  Eq.  (25)  for  change  of  variables,  one  can  finally 
write  down  the  perturbed  beam  current  below 
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6Jb(r)  =  (l/2nbmeSc)  /  dR  — ■  6JQ(r,R) 

“'r 


If  the  plasma  is  a  perfect  conductor  with  infinite  conductivity,  all 
the  beam  particles  which  initially  lie  on  a  magnetic  field  line  should  continue 
to  do  so,  and  when  they  move  the  field  must  in  some  sense  be  carried  along  by 
the  conductor.  This  is  what  Alfven13  expressed  dramatically  by  saying  that 
the  field  and  fluid  were  "frozen  together".  On  the  contrary,  if  the  plasma 
resistivity  is  significant,  there  is  a  phase  lag  between  the  beam  oscillations 
and  the  magnetic  restoring  force.  This  magnetic  decay  force  prevents  the 
field  line  from  following  the  beam  oscillation  which  causes  the  instability  to 
grow.  This  driving  force  is  well  represented  by  the  second  term  of  Eq.  (23). 

The  plasma  return  current  influences  the  structure  of  the  equilibrium  beam 
current  and  vise  versa  which  will  be  discussed  in  a  separate  paper  later 
on.  What  we  concern  here  is  both  currents  with  different  radial  profile 
strongly  effect  the  perturbed  beam  current  which  in  turn  goes  to  Eq.  (23) 
to  determine  the  perturbed  field.  This  is  a  major  driving  force  for  the 
axisymmetric  perturbations.  We  will  concentrate  to  derive  the  most  important 
quantity,  namely  the  perturbed  beam  current  for  the  rounded  beam. 

Recall  that  for  the  loss-cone  distribution  as  in  Eq.  (16)  the  corresponding 
beam  density  profile  is  known  to  be  flat-topped  as  shown  in  Eq.  (17).  To 
solve  the  linearized  Vlasov  equation  for  the  perturbed  distribution  function 
6fb(x,£,t),  one  can  do  the  orbit  integration  along  the  unperturbed  particle 
orbits,  using  the  method  of  characteristic.  From  this  standard  Vlasov  analysis 
and  using  a  variational  approximation  method,  an  approximate  dispersion  relation 
can  be  obtained  with  the  help  of  the  trial  eigenfunction  of  the  form 


n  2  j 

6<t>(r)  «  £  Vr/Rb3 

j-0 
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where  n  is  the  radial  mode  number.  We  will  mainly  be  interested  in  the  lower 
order  number  modes,  namely  n  =  1  and  2  which  are  believed  to  be  the  most  important 
perturbations.  For  n  =  1  which  is  identified  as  sausage  mode,  the  eigenfunction 
has  the  simple  form 


d^jCr)  = 


V1  ■  r2/Rb> 

0 


0  -  r  <  Rb 
*b-  r 


(29) 


which  is  a  self-consistent  solution  to  Eq.  (23).  The  perturbed  beam  current 
density  as  defined  in  Eq.  (24)  can  be  integrated  with  the  help  of  eigenfunction 
(29).  This  standard  procedure  but  tedious  operation  gives  the  final  result 
for  the  axial  component  of  the  perturbed  current  density  like 


5Jb(r)  * 


cel  i  3_  /  Mpb(r) 

9  3r  V 'nTTZ 


3_ 

9r 


(30) 


as  given  in  Eq.  (32)  of  reference  5,  where  the  plasma  frequency  0)^^ 
defined  for  the  flat-topped  beam  by 


(r)  is 


2 

u  ,  ,  ,  * 
pb(r) 


2 

(0)  e  /Ybm 

0 


°<  r  <Rb 

Rb<  r 


(31) 


Similarly,  let  us  consider  next  the  radial  mode  number  n  ■  2  which  is  known 
as  hollowing  mode.  The  name  "hollowing"  comes  because  the  perturbation  can 
indeed  reduce  the  beam  density  on  axis,  leading  to  a  hollow  beam  profile.  The 
eigenfunction  for  n  »  2  has  been  given  by 
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<S4>2  Cr)  = 


A0(l  -  4r2/R^  +  3r4/Rj) 
0 


°<r<Rb 
Rb<  r 


(32) 


After  some  tedious  algebraic  manipulations  as  demonstrated  in  reference  (5) , 
the  z  component  of  the  perturbed  current  density  finally  becomes 


5JbCr) 


-  cr 


4ir 


3 

3r 


(  Ml 

V  3r  Q2  - 


(i  + 


18  4  \ 

n2  -  4u?  7 


(33) 


2  2 

Eqs.  (30)  and  (33)  display  explicitly  the  intrinsic  singularity  at  ft  =  Wg 

for  wave-particle  resonance.  Moreover,  it  is  further  noticed  that  the  hollowing 

mode  has  another  singularity  at  second  harmonic  betatron  frequency  too;  i.e., 

2  2 

ft  *  4o)g.  For  the  same  token,  one  can  expect  more  singularities  at  higher 
harmonic  frequencies  for  higher  order  modes.  Since  u)g  is  single-valued  for  the 
flat-topped  beam,  every  electron  can  be  in  resonance  at  the  same  frequency  and 
there  is  no  limitation  for  the  instability  to  grow.  Therefore,  it  is  almost 
impossible  for  a  flat-topped  beam  to  propagate  in  any  medium.  Fortunately  as  we 
shall  see  later  because  the  beam  is  rounded,  the  nature  is  not  going  to  be  so 
cruel  to  allow  the  instability  to  grow  indefinitely. 

For  the  obvious  reason,  we  are  going  to  consider  the  importance  of  the 
influence  of  rounded  radial  density  profile  on  axisymmetric  perturbations.  As 
we  mentioned  earlier,  the  betatron  frequency  is  a  function  of  radius  r  as  long 
as  the  beam  profile  is  rounded.  It  is  very  common  to  have  energy  spread  for  a 
high  energy  beam  even  right  after  coming  out  of  the  accelerator.  In  nature, 
it  always  represents  an  oscillation  in  energy  about  the  equilibrium  value. 

As  long  as  the  energy  exceeds  the  equilibrium  energy,  the  electron  will  travel 
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in  orbit  of  larger  radii,  conversely,  it  will  have  smaller  orbit  at  lower  energy. 
The  electrons  reside  on  each  infinitesmal  Ar  can  have  approximately  constant  u)g 
with  which  they  resonate  locally.  Therefore,  it  is  quite  reasonable  to  extend 
from  the  previous  results  for  the  flat-topped  beam  profile  to  include  the 
rounded  beam  profile  effect.  To  investigate  the  sausage  and  hollowing  modes  for 
the  rounded  beam,  we  have  to  derive  first  the  perturbed  beam  current  and  the 
associated  eigenfunction  for  each  mode  which  are  described  below  separately. 

1.  Sausage  mode  (n  =  1) 

It  is  evident  that  a  beam  with  rounded  radial  density  profile  can 
be  subdivided  into  groups  of  particles  with  the  same  value  of  energy,  and  in 
addition,  it  can  also  be  sliced  into  components  with  different  radii.  As  pointed 
out  in  Eq.  (30) ,  each  energy  group  designated  by  R  gives  its  perturbed  axial 
current  density  for  the  n  =  1  sausage  mode 


6J0(r,R) 


27re2cS 


2  1  a_ 

t  3r 


dJ,  (R)/dR 

D 

n2  -  u>2(r) 


(34) 


where  5JQ(r,R)  should  be  a  continuous  function  on  0  (  r  (  R  and  vanish  smoothly 
at  r  =  0.  Substituting  Eq.  (34)  into  Eq.  (27)  for  summing  up  each  individual 
component  of  the  perturbed  beam  current,  we  have  finally  the  result  for  the 
axial  perturbed  beam  current  like 


,T  ,  A  eg  1  9 

oJ,  (r)  *  —  —  r — 

b  ym  r  5r 


dJ.  (R)/dR 


/  dR 

JT  ST  -  w 


_  3_£i> 
2  9r 

8 


The  corresponding  eigenfunction  is  given  by 
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r* 


V 


Note  that  the  intrinsic  singularity  at  ft  =  wg  in  Eq.  (30)  for  the  flat- topped 
case  has  been  replaced  by  the  integration  over  R  space  which  eventually  gives 
a  logarithmic  singularity.  Because  of  this,  the  growth  rate  is  no  longer 
unbounded  as  in  the  flat-topped  case.  The  natural  oscillation  frequency  wg 
varies  for  different  electrons  at  different  radius  r,  so  that  there  is  no 
single  frequency  0  for  which  the  entire  beam  electrons  can  resonate  with 
the  wave. 

2.  Hollowing  mode  (n  =  2) 

Accordingly,  after  considering  the  n  =  1  sausage  mode,  likewise 
we  can  easily  derive,  based  on  the  same  theory,  the  perturbed  current  density 
and  corresponding  trial  eigenfunction  for  the  n  =  2  hollowing  mode.  By  use 
of  Eqs .  (32)  and  (33),  one  can  obtain  the  perturbed  beam  current  which  has  the 
axial  component 


6Jb(r) 


ym  r  3  r  (  Jt 


dJ  (R)/dR 

dR  —2 - —  (1  + 

2  2 
n  -  uig 


18  u 


8 


a2-  4 


■>  m 


(37) 


and  the  corresponding  eigenfunction  for  n  =  2  becomes 


6<t>2(r)  dR  aq(1  -  4r2/R2  +  3r4/R4)  dJb(R)/dR 


(38) 


Obviously,  if  the  beam  current  Jb(R)  is  a  step  function,  then  the  above 
derived  results  for  the  rounded  beam  go  back  to  that  for  the  flat-topped 
beam  case. 

B.  Dispersion  relation 

The  derived  eigenvalue  equation  (23)  along  with  the  perturbed  beam 
current  (35)  or  (37)  constitute  the  main  result  for  axisymmetric  perturbat.  .'ns 
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which  can  be  summarized  in  the  general  form 


r  Ir  [r  <Ir  5l^r)  Fj  +  °(r)wi6$(r)  =  0 


where  j  =  1  or  2  representing  sausage  and  hollowing  mode  respectively,  so  that 
the  function  I\  (r)  are  defined  by 


r,(r)  .  1  „  iHSi  fiSL 
1  ymc  J 


4-rreS  f"  dJbW/dR 


ymc  J  2  2 

*r  0  -  u6(R) 


4-neB  f 

JT 


dJu(R)/dR  18u‘(R) 

dR  ~2  2  ~~2  2 
ft  -  u>g(R)  ST  -  4ujgCR) 


It  is  worth  while  to  point  out  that  the  eigenvalue  problem  represented  by 
Eq.  (39)  is  independent  of  o(0),  J^CO)  and  3  which  enables  us  to  use  them  as 
a  normalization  factor  for  w  and  ft  later  in  expressing  the  dispersion 
relationship.  Therefore,  the  dispersion  relation  can  be  expressed  quite 
generally  as  a  functional  relationship  between  ft  /u^B  and  where  the 

•7  7 

magnetic  decay  time  for  the  perturbed  current  is  defined  by  =  ira(0)R^/2c  . 
For  any  rounded  beam,  Eq.  (39)  along  with  (40)  or  (41)  can  be  used  to 
investigate  the  influence  of  the  magnetic  decay  time  and  the  plasma  return 
current  on  axisymmetric  instabilities. 

Ironically,  Eq.  (39)  is  a  standard  eigenvalue  problem  which  can  be 
solved  numerically.  However,  the  method  we  adapted  here  enables  us  to 
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study  certain  specific  cases  which  exhibit  explicitly  more  physics.  The 
dispersion  relation  that  determines  the  complex  eigenfrequency  w  can  be 
obtained  easily  by  multiplying  Eq.  (39)  by  r  (r)  on  both  sides  and 
integrating  over  r  from  0  to  «  .  This  operation  gives  the  following  result 


oo  oo 

Jq  dr  r  j^d5<j)(r)/drj2  I\(r)  =  j^dr  r  a(r)j^<54>(r)J 


(42) 


The  dependence  of  stability  properties  on  the  form  of  the  conductivity 
profile  is  of  considerable  interest.  Since  the  right  hand  side  of  Eq.  (42) 
is  not  a  function  of  0,  the  electrical  conductivity  a  (r)  only  serves  as  a 
weighted  function  in  determining  the  value  of  u.  In  other  words,  the  results 
become  free  of  restriction  to  beams  propagating  in  a  medium  with  any  fixed 
conductivity  profile.  For  any  growing  wave,  broadening  the  conductivity 
channel  (R^  >  R^)  will  reduce  u  and  has  a  stabilizing  effect.  A  narrow 
profile  for  o(r)  (R^<  ^  can  substantially  enhance  the  growth  rate. 
However,  the  equilibrium  conductivity  c(r)  in  reality  is  a  response  to  beam 
current  and  ambient  plasma  according  to  experimental  condition  and  must 
be  allowed  to  change  with  respect  to  time  and  space.  We  will  examine  this 
important  property  in  the  near  future;  Eq.  (42)  can  be  solved  numerically 
to  determine  the  stability  properties  for  various  system  parameters  of  our 
interest.  However,  if  the  beam  current  is  a  step  function  as  in  the 
flat-topped  beam  case,  Eq.  (42)  could  lead  to  the  following  analytical 
results.  Substituting  Eqs.  (36)  and  (40)  into  Eq.  (42),  the  dispersion 
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relation  for  the  m  =  1  (sausage  mode)  becomes  the  simple  form 


4  . 

3  luTd  ■  ,2 


"„b8t  »  -  2f> 

2  >6b  »  *  f> 


Similarly,  by  use  of  Eqs.  (38),  (41)  and  (42),  the  dispersion  relation  for  the 
n  =  2  (hollowing  mode)  becomes 


202 
u>  8,, 


?  2 
9  u  ref 


8  -2>Bb  ^  L  8  -  8  vBb  (l  - f) 


Since  the  right  hand  side  of  Eqs.  (43)  and  (44)  are  real  always,  we  can  expect 
Re  (w)  =  0  all  the  time.  So  the  instabilities  for  the  flat-topped  beam  are 
purely  growing  modes. 

C.  Numerical  results 

In  the  dispersion  relations  (42) ,  we  fix  the  real  frequency  12  and 

solve  for  the  complex  u  numerically  for  various  system  parameters  to  determine 

the  wave  properties.  The  radial  dependence  of  the  equilibrium  beam  current, 

steady  plasma  current  and  conductivity  of  plasma  channel  are  all  assumed  to 

have  the  similar  Bennett  profile,  but  different  radii  are  allowed  in  general 

to  express  the  complexity  of  the  nature.  In  the  following,  the  sausage  and 

hollowing  modes  will  be  investigated  separately  for  various  system  parameters 

such  as  f,  and  R^.  The  dispersion  relation  will  be  expressed  as  a 

2  2  2 

functional  relationship  between  fi  and  The  results  between 

rounded  and  flat-topped  beams  will  be  compared  with  physical  interpretation. 
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1.  Sausage  mode  (n  =  1) 

Let  us  start  considering  the  n  =  1  sausage  mode  in  the  simplest 

case  where  =  R^;  i.e.,  the  beam  and  plasma  return  currents  have  the  same 

Bennett  profiles  of  equal  radii  so  that  f  =  Jp(r)/J^(r)  is  constant  everywhere 

for  any  radius  r.  According  to  Eq.  (36),  if  the  beam  current  J,  (R)  has  the 

b 

Bennett  profile  as  expressed  in  Eq.  (1),  then  64^ (r)  can  be  integrated  easily 
to  obtain 


fi*1(r)  =  (1  +  R^/r2)/2  -  ln(l  ♦  R^/r2)-0.S/  (1+R^/r2) 

Rb 


(45) 


which  has  the  peak  value  at  r  =  0  and  decreases  monotonically  *0  zero  at  r  =  » 
Furthermore,  the  eigenfunction  C^^(r)  and  its  first  derivative  with  respect  to 
r  (i.e.,  3^/3^  are  independent  of  ft  and  f,  so  that  the  LHS  of  dispersion 
relation  (42)  decides  the  stability  criteria  (stable  or  unstable)  and  the  RHS 
determines  the  strength  of  the  instability.  Using  Eqs.  (1)  and  (11),  Eq.  (.40) 
can  be  integrated  to  give  the  following  analytical  form 


r 

‘  1 


(r) 


_j _ +  n2 

1.4'  >6b2“-f) 

"b 


lnfl 


(1^)  a 

Rb 


2  J 


(i-f) 


-1 


(46) 


Note  that  ^(r)  is  complex  with  the  definition  for  the  complex  logarithm 

In  2  i  In  r  ♦  i6  (47) 

2  2  1/2  -1 

where  z  =  x  ♦  iy  and  r  =  (x  +  y  )  ,  0  =  tan  (y/x) .  It  becomes  apparent 

that  the  intrinsic  singularity  appears  in  Eq.  (30)  has  been  replaced  by 
the  logarithmic  singularity  as  explicitly  shown  in  Eq.  (46)  as  we  mentioned 
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before;  for  any  fixed  ft  and  f,  the  pole  for  logarithmic  singularity  can  be 

located  from  the  limit  lim  £n(z)  which  gives 

z-+0 


It  is  because  5<j>j (r)  obtained  in  (45)  is  a  strong  monotone  decreasing  function, 
with  respect  to  increasing  r,  the  integrant  in  the  LHS  of  Eq.  (42)  is 
almost  zero  near  the  logarithmic  singularity  Rg.  In  this  regard,  the 
perturbation  is  contributed  essentially  from  the  integration  interval 
between  0  and  Rg.  Otherwise  the  LHS  of  Eq.  (42)  could  approach  infinity 
which  gives  the  unbounded  growth  rate. 

The  results  of  solving  Eqs.  (42)  and  (45)  give  the  unstable  regions  as 
shown  in  Fig.  1  for  different  values  of  f  and  ft.  The  dashed  and  solid  lines 
are  for  the  flat-topped  and  rounded  beams  respectively.  As  Eq.  (43) 
indicated,  the  unstable  region  for  the  flat-topped  beam  is  bounded  by  two 
straight  lines  in  the  f  -  ft  space.  The  value  of  f  is  restricted  but  taken 
from  -1  to  1  which  covers  completely  the  possible  magnitude  and  direction  of 
the  plasma  return  current  in  our  case.  The  fact  is  for  a  single  pulse 
electron  beam,  the  plasma  return  current  is  known  to  be  flowing  in  the 
direction  opposite  to  the  beam  current.  However,  the  situation  is  not  always 
so  simple.  For  example,  as  in  the  case  where  the  beam  is  a  continuous  pulse 
train,  the  intermission  between  the  pulse  can  cause  the  plasma  return  current 
going  in  the  opposite  direction  before  the  second  pulse  arrives.  Anyway,  the 

unstable  region  for  the  rounded  beam  has  been  greatly  reduced  towards  smaller 
n  in  the  upper  half  of  f  space.  It  is  interesting  to  observe  that  the 
instability  exists  even  when  ft  =  0  for  both  beams.  For  ft  =  0,  Eq.  (46) 
can  be  simplified  to  give 
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rx  (r)  =  1  -  (1  ♦  r2/!^)"1  (1  -  f)'1  (49) 

The  real  frequency  Re(w)  =  0  for  the  rounded  beam  also  indicating  that  the 
perturbations  are  purely  growing  modes,  which  can  only  be  driven  unstable  by 
certain  amount  of  plasma  return  current,  i.e.,  f  ^  0.  As  f  increases,  the 
rounded  beam  becomes  unstable  first  at  f  *  0.2  approximately  while  the  flat- 
topped  beam  remains  stable  until  about  f  *  0.5.  The  growth  rates  Im(wr^)  for 
both  beams  are  plotted  against  f  in  Figure  2  when  Q  =  0.  The  growth  rates  are 
approaching  infinity  in  the  limit  of  f  =  1  as  one  can  expect.  Because  the 
plasma  return  current  cancels  the  beam  current  completely  and  the  self-magnetic 
pinched  force  no  longer  exists. 

2 

If  fl  ^  0,  the  growth  rates  versus  0  for  the  flat-topped  beam  are  shown 

in  Figure  3  for  various  values  of  fractional  current  neutralization  f.  The 

growth  rates  could  go  to  infinity  no  matter  what  value  of  f  is  chosen  because 

of  the  singularities  which  appear  explicitly  in  Eq.  (43).  Again,  the  results  of 

Re(oj)  are  always  zero  indicating  that  the  perturbation  are  purely  growing  modes. 

As  to  the  results  for  rounded  beam,  the  ImCux^)  and  Re(cord)  for  various  f  are 
2 

shown  versus  Q  in  Figure  4(a)  and  (b)  respectively.  The  growth  rate  is  very 
small  when  f  =  0  where  the  sausage  mode  is  driven  by  magnetic  decay  force  alone. 

For  larger  f,  the  growth  rate  increases  dramatically  which  indicates  the  strong 
influence  of  the  plasma  return  current  on  the  instability.  In  comparison,  the 
rounded  beam  does  not  have  the  feature  of  unbounded  linear  growth  rate  as  in 
the  case  for  the  flat-topped  beam.  This  is  an  important  quantity  for  the  rounded 
beam  which  is  in  favor  of  the  beam  propagation.  As  we  know,  if  the  beam  is  rounded, 
then  the  betatron  frequency  for  each  electron  is  not  unique.  A  wave  with  a  given 
real  frequency  (2  is  resonant  with  portion  of  the  beam  electrons  only,  but  not 


h 
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all  of  them,  and  the  resonant  electrons  tend  to  be  localized  in  particular 

radial  ranges  of  the  beam;  the  oscillation  frequency  Re(u)  as  shown  in 

2 

Figure  4(b)  increases  linearly  with  respect  to  n  .  The  reason  of  being 

that  is  very  simple.  As  one  can  see  in  Eq.  (42)  the  Re(w)  comes  from 

Imd^r)),  and  the  imaginary  part  of  T^(r)  is  contributed  only  from  the  pole 

at  0,  ~  ujg  according  to  Eq.  (40)  if  the  path  of  integration  is  chosen  along  the 

Re(R)  axis.  By  Cauchy's  theorem  a  contour  can  be  deformed  without  crossing  the 

singularity  and  then  applying  the  theorem  of  residues  for  evaluating  the 

definite  integral  in  Eq.  (40),  which  gives  eventually  the  same  unique  value  of 

iri.  More  specifically,  the  last  term  in  Eq.  (46)  has  a  logarithmic  singularity 

from  which  a  value  of  iri  is  obtained  according  to  Eq.  (47).  Clearly,  the  term 

2 

is  linearly  proportional  to  ft  . 

The  electrical  conductivity  of  the  plasma  channel  is  a  rather  complicated 
quantity  characterized  by  a  radial  profile  o(r)  which  determines  the  radial 
distribution  of  the  plasma  return  current.  For  simplicity,  the  electrical 
conductivity  o(r)  and  the  plasma  return  current  J^(r)  are  assumed  to  be  steady 
and  have  the  Bennett  profiles  of  equal  radii;  usually  the  plasma  return 
current  flows  at  larger  radius  than  the  beam  current  at  lower  pressure  because 
electric  breakdown  considerably  spreads  the  conductivity  channel.  On  the 
contrary,  the  higher  pressure  forces  the  plasma  return  current  concentrating 
towards  the  beam  axis.  Therefore,  without  loss  of  generality,  it  is  our  interest 
to  consider  these  two  extreme  cases  where  R^  j*  R^,  i.e.,  the  plasma  return 
current  J  (r)  and  beam  current  J,  (r)  have  the  same  Bennett  profiles  but 

p  D 

different  radii.  When  we  do  this,  it  is  important  to  remember  that  the 
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restriction  on  fractional  current  neutralization  f  as  shown  in  Eq.  (13)  has  to 

be  .enforced  so  that  the  equilibrium  beam  profile  may  exist.  Furthermore, 

under  the  condition  of  Rp  ^  R^,  it  is  more  convenient  to  use  the  effective 

current  neutralization  fg££  as  defined  in  Eq.  (5)  instead  of  f. 

First,  let  us  consider  the  case  of  Rp  >  R^.  Figure  (5)  demonstrates 

the  change  of  unstable  region  in  the  (1  -  fg£^  space  due  to  different  ratio  of 

Rp/Rb.  As  feff  =  0,  the  unstable  regions  in  £2  domain  are  identical  because 

there  is  no  plasma  return  current  so  that  the  ratio  of  Rp/Rb  is  irrelevant. 

The  electrical  conductivity  a(r)  serves  as  a  weighting  function  in  determining 

the  growth  rate  of  instability.  For  fg££  $  0,  the  instability  growth  region 

shrinks  in  the  £2  space  but  extends  towards  higher  negative  fg££  region  as  the 

ratio  of  Rp/Rb  increases.  If  fg££  <  0,  the  expanded  plasma  return  current 

flows  in  the  same  direction  as  the  beam  current  does,  which  is  contrary  to 

the  concept  used  in  the  center-core  transmission  line  where  both  currents  flow 

in  the  directions  opposite  to  each  other.  No  wonder  it  is  more  unstable  for 

larger  ratio  of  Rp/Rb;  as  far  as  the  growth  rate  is  concerned,  Figure  (6)  shows 

2 

the  solution  (wr^)  versus  £2  for  various  fg££  when  Rp  =  2  R^.  The  growth  rates 

ImCojT^)  increases  as  fg££  increases.  Obviously,  the  growth  rates  and  the 

magnetic  decay  force  ReCcox^)  are  smaller  when  comparing  to  those  in  Figure  (4) 

where  R  =  R,  ,  which  indicates  that  the  increase  of  plasma  current  radius  aids 
p  o 

stability  for  feff  >  l’ 

Now  let  us  consider  the  other  extreme  case  where  Rp  <  Rb,  the  results  for 
the  unstable  regions  in  £2  -  fg££  space  are  shown  in  Figure  (7)  for  different 
ratio  of  Rp/Rb.  The  current  neutralization  f  or  effective  current 
neutralization  fg££  has  to  be  restricted  according  to  Eq.  (14)  as  a 
minimum  requirement  for  the  radial  confinement  of  the  Bennett 
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beam.  As  a  result,  there  is  a  forbidden  region  for  each  case  in  Figure  (7) 

2 

which  fulfills  the  inequality  fg££  <  (Rp/R^)  in  fg£^  domain;  e.g., 
fgff  >  .25  when  Rp  =  0.5  R^. 

As  the  plasma  channel  is  narrower  than  the  beam  current,  the  unstable 
region  increases  in  fl  space  for  fg^£  >  0  but  remains  very  stable  for  fgf£  <  0. 
As  fg££  <  0,  the  self-magnetic  pinched  force  has  been  enhanced  somewhat  by  the 
plasma  return  current  which  flows  in  the  same  direction  as  the  beam  current 

so  that  the  beam  is  very  stable.  The  growth  rates  for  Rp  =  0.5  Rb  are  shown 

~> 

in  Figure  (8)  versus  (T  for  different  fg££.  The  value  of  fg££  is  chosen  only 
from  0  to  0.2  which  corresponds  that  f  increases  from  0  to  0.8  because  of  the 
restriction  of  Eq.  (14).  The  growth  rate  can  be  very  large  even  for  small 
value  of  fg££.  Finally,  let  us  summarize  the  sausage  mode  calculations  by 
presenting  the  criteria  for  non-oscillatory  growth  region  in  the  fgf£  -  Rp/Rb 
space  which  is  shown  in  Figure  (9) .  It  is  a  purely  growing  mode  where  (2=0 
and  Re  (to)  =  0.  As  one  can  see,  the  sausage  mode  of  this  kind  can  only  be 
excited  when  fg££  >  0. 

2.  Hollowing  mode 

Next  to  the  sausage  mode,  we  want  to  consider  and  therefore  do 
the  similar  calculations  for  the  hollowing  mode  (n  =  2).  For  the  simplest  case 
where  Rp  =  Rb,  the  eigenfunction  (38)  has  the  form  after  integration 


5<f>2(r)  =  64>-, (r)  -  1 . 5(r/ F^) 4  |o. 75  ♦  0.25 


X  -  1.5  X 


+  1 . 5  1  n  X  + 


(50) 
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where  ^(r)  is  equivalent  to  that  shown  in  Eq.  (45)  and  X  =  1  +  R^/r2  for 
simplification.  6$2(r)  is  also  a  monotone  decreasing  function  which  vanishes 
as  r  -*■  ».  Likewise  the  expression  for  r2(r)  in  Eq.  (41)  can  be  integrated 
to  give 


r2(r)  -  Fj  (r)  ♦  -J- 
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(51) 


where  ^(r)  denotes  the  same  expression  as  in  Eq.  (46)  and  X  =  1  +  (R^r)2 

also  for  simplification.  As  one  can  see  in  Eq.  (51),  the  last  two  terms 

contain  two  logarithmic  singularities  which  have  been  formed  due  to  the  two 

intrinsic  singularities  in  Eq.  (41).  Substituting  Eqs.  (50)  and  (51)  into 

Eq-  (-*2),  we  solve  the  dispersion  relation  numerically  for  various  system 

parameters.  The  results  for  the  unstable  regions  in  the  £1  -  f  space  are 

illustrated  in  Fig.  (10)  for  the  flat-topped  (dashed  line)  and  rounded  beam 

(solid  line)  respectively.  Note  that  the  unstable  region  for  the  flat-topped 

beam  is  greatly  reduced  by  the  influence  of  the  rounded  beam  profile.  From 

Eq.  (44)  the  growth  rates  for  the  flat-topped  beam  can  be  obtained  easily 

and  shown  in  Fig.  (11)  for  various  value  of  current  neutralization  f. 

Correspondingly,  the  Im(u)T^)  and  ReCwr^)  for  the  rounded  beam  are  plotted 
2 

against  ft  for  various  f  in  Fig.  (12).  As  one  can  immediately  see,  the 

unbounded  linear  growth  rates  of  the  flat-topped  beam  are  replaced  by  the 

finite  growth  rates  of  the  rounded  beam.  For  the  rounded  beam,  both  Im(u)T.) 

d 

and  Re(uTcj)  increases  as  f  increases.  Note  that  the  growth  rate  is  very 
large  even  when  f  =  0,  the  hollowing  mode  is  driven  unstable  by  the  magnetic 
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decay  force  along.  It  is  further  noticed  that  the  peaks  of  the  growth  rates 
can  be  found  near  the  dashed  Kne  as  shown  in  Fig.  (10)  where  the  growth 
rates  for  the  flat-topped  beam  become  infinity. 

Now,  we  are  ready  to  make  a  comparison  between  the  sausage  and  hollowing 
mode  perturbations  for  the  rounded  beam.  The  growth  rates  as  shown  in  Fig.  (12) 
for  the  hollowing  mode  are  generally  larger  than  that  of  the  sausage  mode 
in  Fig.  (4).  Fow  example,  when  f  *  0,  the  growth  rate  for  the  hollowing  mode 
is  at  least  an  order  of  2  larger  than  the  n  =  1  sausage  mode.  It  concludes 
that  the  magnetic  decay  force  plays  an  important  role  in  exciting  the  hollowing 
mode.  As  f  increases,  the  ratio  of  the  growth  rate  between  the  two  modes 
reduces  gradually  but  still  remains  a  factor  of  10  approximately  when  f  =  0.8. 

As  f  <0,  the  plasma  return  current  has  the  stabilising  effect  on  both  modes, 
because  it  introduces  additional  self-magnetic  pinched  force  for  the  beam. 

The  sausage  modes  are  very  stable  but  not  for  the  hollowing  modes  because  the 
growth  rates  are  very  large  even  when  f  =  0,  it  certainly  needs  more  plasma 
return  current  flowing  in  the  same  direction  to  stabilise  the  beam.  Generally 
speaking,  because  of  the  huge  growth  rate  of  the  hollowing  mode,  it  seems 
unlikely  to  have  any  sausage  mode  oscillations  or  enough  time  for  them  to  grow 
before  the  hollowing  mode  totally  destroys  the  beam.  Now,  if  one  looks  closely 
in  Figs.  (4)  and  (12)  again,  one  can  find  another  big  difference  between  the 
results  for  the  sausage  and  hollowing  modes.  Namely  the  maximum  growth  rate 
for  various  f  occurs  near  small  value  of  fi  for  the  sausage  mode  but  span  to 
very  large  values  of  £1  except  for  those  with  high  current  neutralication 
for  the  hollowing  mode.  Recall  that  the  perturbation  we  are  interested  is  in 
the  low  frequency  regime.  The  theory  derived  here  is  not  necessarily  correct 
for  high  frequency  perturbation.  The  beam  electrons  at  various  energy  level 
oscillate  with  various  betatron  frequency  at  different  radius.  From  Cq.  (11), 
the  betatron  frequency  u)g(r)  has  the  maximum  value  u^(0)  at  the  beam  axis  and 
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decreases  outwards  outwards  radially.  Because  of  the  natural  oscillation  of  the 
electrons,  it  is  likely  that  the  waves  are  excited  primarily  in  the  frequency 
range  of  0  <  ft  <  w  (0) .  In  this  respect,  although  the  growth  rates  of  the 

P 

hollowing  mode  are  very  large  in  general,  one  has  to  be  very  careful  about 
th-e  location  of  Q  where  the  peaks  of  growth  rates  occur.  One  could  argue 
that  it  is  almost  impossible  to  have  waves  excited  with  very  large  frequency 
ft  (e.g.  0  >  cog  CO)  )  unless  the  perturbation  is  imposed  externally.  In  Fig.  (12) 
the  maximum  growth  rate  occurs  near  the  region  where  ft  >  Wg(0)  which  is  not 
likely  to  happen  in  the  nature.  Nevertheless,  the  maximum  growth  rate  for 
f  >  0.5  can  be  found  in  the  region  of  ft  satisfying  ft  <  wQ(0).  It  indicates 

P 

the  possible  excitation  of  the  hollowing  mode  naturally  when  the  current 
neutralization  is  high.  There  has  been  some  evidence  from  numerical  simu- 

g 

lation  results.  For  example,  there  is  no  sausage  mode  oscillation  which 
has  even  been  observed  experimentally  for  a  self-pinched  electron  beam 
propagating  in  the  plasma.  In  addition,  numerical  simulation  has  indicated 
that  the  electron  beam  becomes  hollow  if  the  current  neutralization  exceeds 
0.5  or  so. 

It  is  straight  forward  to  extend  the  calculation  of  hollowing  mode  to 

consider  the  case  where  f  R^,  namely  the  beam  and  plasma  return  currents 

have  the  similar  Bennett  profiles  but  unequal  radii.  First,  for  the  case 

Rp  >  Rb,  Fig.  (13)  illustrates  the  various  growth  regions  in  the  ft  -  feff 

space  for  different  ratio  of  Rp/R^.  The  beam  seems  more  stable  as  the 

ratio  Rp/Rfe  increases,  because  the  unstable  region  in  ft  space  have  been 

greatly  reduced  and  concentrated  towards  higher  ft  domain.  The  growth  rates 
2 

versus  ft  for  Rp  =  2R^  are  shown  in  Fig.  (14)  for  various  fe^.  One  can  see 
clearly  that  the  growth  rates  are  not  a  strong  function  of  f  ££  any  more. 

The  reason  is  very  simple,  the  hollowing  mode  perturbation  is  highly  unstable 
at  f  *  0  where  the  instability  is  driven  simply  by  the  magnetic  decay  force. 
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The  plasma  return  current  tries  to  change  the  situation.  However,  under  the 

condition  of  R^  >  R^,  the  amount  of  plasma  return  current  overlapped 

with  the  beam  remains  almost  unchanged  even  for  very  large  ratio  of  R  /R^. 

On  the  other  hand,  for  the  second  case  where  R^  <  R^,  the  unstable  region 
for  the  hollowing  mode  perturbation  as  shown  in  Fig.  (15)  expands  very  much 
in  2  space  as  the  ratio  of  R^/R^  decreases  first  (e.g.,  R^  =  0.5  R^) ,  but 
shrinks  back  for  further  decrease  of  R^/R^,  a  forbidden  region  in  the  f^^ 
domain  is  drawn  in  Fig.  (15)  where  the  equilibrium  Bennett  beam  fails  to 
exist.  The  growth  rates  for  R^  =  0.5  R^  are  plotted  in  Fig.  (16)  against 
fs”  for  various  f  It  is  very  clear  that  as  f  ^  increases  the  growth 

rate  increases  and  car.  go  to  infinity  (e.g.,  f  ^  >  0.15  or  f  >  0.6 
equivalently).  For  higher  fg££.  this  is  a  very  unstable  situation  for  the 
beam  because  of  the  hollowing  mode  oscillations.  Once  the  self-magnetic 
pinched  force  is  reduced  near  the  axis  by  the  presence  of  the  plasma  return 
current,  it  is  not  difficult  to  understand  why  the  beam  electrons  on  the 
beam  axis  can  easily  be  expelled  by  the  highly  concentrating  return  current. 
As  long  as  the  hollowing  mode  exists,  the  beam  is  going  to  become  annular 
eventually  and  expand  continuously. 
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IV.  Conclusions 

The  equilibrium  and  stability  of  an  intense  relativistic  self-pinched 
electron  beam  propagating  in  a  resistive  plasma  channel  has  been  studied  by 
the  Vlasov -Maxwell  Theory.  For  the  thermal  equilibrium  distribution  function, 
the  steady  and  equilibrium  electron  beam  without  plasma  return  current  has 
the  familiar  Bennett  radial  density  profile.  When  the  plasma  return  current 
is  unavoidable  in  the  system,  the  minimum  requirement  for  the  radial  confinement 
of  the  electron  beam  has  been  derived  to  restrict  the  degree  of  current 
neutralization.  An  analytical  theory  of  the  axisymmetric  perturbations  including 
the  sausage  (n  =  1)  and  hollowing  (n  =  2)  modes  has  been  derived  for  an  electron 
beam  with  rounded  radial  current  profile.  Complete  space-charge  neutralization 
and  the  paraxial  beam  electron  motion  have  been  assumed  in  the  derivation. 

In  this  paper,  we  have  considered  in  particular  a  wave  with  real  frequency 
and  Im(oj)  >0  which  has  a  certain  wave  amplitude  at  the  beam  head  and  grows 
spatially  backward  along  the  beam.  It  is  essentially  a  convective  instability, 

A  numerical  calculation  has  been  made  for  the  special  case  where  the  equilibrium 
radial  distribution  of  the  beam  and  plasma  return  current,  as  well  as 
electrical  conductivity  of  the  plasma  channel,  have  similar  Bennett  profiles 
but  different  radii  is  also  allowed.  The  results  for  axisymmetric  perturbations 
have  shown  that  the  important  system  parameters  such  as  the  effective  current 
neutralization,  the  radial  profile  of  electrical  conductivity  of  plasma  channel, 
beam  and  plasma  return  current  are  very  sensitive  in  determining  the  stability 
behavior  of  an  electron  beam.  First,  the  results  for  the  rounded  beam  have  been 
compared  to  that  of  the  flat-topped  beam.  Generally  speaking,  there  is  no  linear 
bound  on  wave  growth  for  the  flat-topped  beam.  Fcr  the  rounded  beam,  the  finite 
growth  rates  for  the  hollowing  mode  are  much  larger  than  that  of  the  sausage 
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mode.  In  practice,  it  has  been  concluded  that  the  naturally  exited  hollowing 
mode  oscillation  is  more  likely  to  happen  when  the  current  neutralization 
exceeds  about  0.5. 

It  has  been  shown  from  the  calculations  the  high-frequency  hollowing  mode 
2  2  2 

perturbations  (i.e.,  ft  can  be  driven  unstable  by  the  magnetic  decay 

force  alone  where  the  phase  lag  between  the  transverse  displacement  of  the 
beam  and  the  magnetic  field  drives  the  instability.  The  naturally  occurring 
plasma  return  current  can  flow  either  in  the  same  or  in  the  opposite  direction 
as  that  of  the  beam  current.  The  one  flows  in  the  same  direction  introduces 
additional  self-magnetic  pinched  force  which  intends  to  stabilize  the  beam. 

On  the  contrary,  the  one  in  the  opposite  direction  reduces  the  pinched  force 
and  thus  enhances  the  perturbations.  As  R^  <  R^  where  a  narrow  conducting 
channel  confines  the  plasma  return  current,  the  highly  concentrating  return 
current  on  beam  axis  is  very  dangerous  for  hollowing  and  destroying  the  beam. 

Finally,  it  is  of  importance  to  compare  the  growth  rates  of  the  axisymmetric 
modes  to  those  of  the  hose  instability.  Generally,  the  growth  rates  for  the 
sausage  mode  are  much  smaller  than  that  of  the  hose  modes  as  obtained  from 
Figure  5  of  Reference  2.  This  may  well  explain  why  the  sausage  mode  oscillation 
has  not  been  found  experimentally  or  in  numerical  simulation.  It  is  simply 
because  of  the  domination  of  the  hose  mode,  the  beam  could  Le  expelled  from 
its  channel  even  before  the  sausage  mode  starts  to  grow.  However,  the  maximum 
growth  rates  for  the  hollowing  mode  are  larger  in  most  cases  or  at  least 
comparable  to  those  of  the  hose  mode  which  addresses  very  clearly  the  importance 
of  the  hollowing  mode  oscillation  in  the  beam  propagation. 
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FIGURE  3.  THE  GROWTH  RATE  OF  THE  SAUSAGE  MODE  VERSUS  ft2  FOR  VARIOUS  f  FOR  THE 
FLAT-TOPPED  BEAM 
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FIGURE  5. 


THE  UNSTABLE  REGION  OF  THE  SAUSAGE  MODE  OF  THE  BENNETT  BEAM  IN  THE 
feff  -  n  SPACE  FOR  Rp/Rb  »  2,  4,  6  RESPECTIVELY 
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FIGURE  7. 
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THE  UNSTABLE  REGION  OF  THE  SAUSAGE  MODE  OF  THE  BENNETT  BEAM  IN  THE  f 
SPACE  FOR  R_/Rb  *  0.9,  0.5,  AND  0.1,  RESPECTIVELY 
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FIGURE  9.  NONOSCILLATOR Y  GROWTH  REGION  FOR  A  BENNETT  BEAM  IN  THE  feff  -  Rp/Rb  SPACE 
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FIGURE  12.  GROWTH  RATE  lm(w)  AND  REAL  EIGENFREQUENCY  RB(w)  VERSUS  P2  FOR  THE  HOLLOWING 
MODE  OF  A  BENNETT  BEAM  WITH  Rp  =  Rb  FOR  VARIOUS  f 
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FIGURE  13.  THE  UNSTABLE  REGION  OF  THE  HOLLOWING  MODE  OF  THE  BENNETT  BEAM  IN  THE 
feff  -  n  SPACE  FOR  Rp/Rb  =  1.1.  2,  4,  RESPECTIVELY 
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FIGURE  14.  GROWTH  RATE  lm(tj)  AND  REAL  EIGENF RECUE NCY  Re(w)  VERSUS  n2  FOR  THE 
HOLLOWING  MODE  OF  A  BENNETT  BEAM  WITH  Rp  =  2  Rb  FOR  VARIOUS  fgff 
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FIGURE  15.  THE  UNSTABLE  REGION  OF  THE  HOLLOWING  MODE  OF  THE  BENNETT  BEAM  IN  THE 
feff  -  fi  SPACE  FOR  Rp/Rb  =  0.9,  0.5,  0.1,  RESPECTIVELY 
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